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is a sum of generators of gauge symmetries, the foliation-preserving diffeomorphisms, and vanishes 
on shell. The scalar constraint is second class, except for a global, first-class part that generates 
time reparametrizations. A reduced phase space formulation is given in which the local part of the 
scalar constraint is solved formally for the lapse as a function of the 3 metric and its conjugate 
momentum. In the infrared limit the scalar constraint is linear in the square root of the lapse. 
For asymptotically flat boundary conditions the Hamiltonian is a sum of bulk constraints plus a 
boundary term that gives the total energy. This energy expression is identical to the one for Einstein- 
aether theory which, for static spherically symmetric solutions, is the usual Arnowitt-Deser-Misner 
energy of general relativity with a rescaled Newton constant. 
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I. INTRODUCTION 

Hofava [1[ has proposed a theory of gravity that 
is closely related to general relativity (GR), but is 
power-counting renormalizable and possesses a pre- 
ferred spacelike foliation of spacetime that breaks the 
spacetime diffeomorphism symmetry down to time- 
dependent 3-dimensional diffeomorphism symmetry and 
a reparametrization of global time. Several variants of 
this theory have been considered (see Ref. @ for a re- 
view). In this paper we study the "consistent extension" 
of Bias, Pujolas and Sibiryakov Q ■ This theory is just the 
"nonprojectable" version of Hofava's original proposal 
with the inclusion of terms involving (lniV)^, the spa- 
tial gradient of the log of the lapse function N, which are 
compatible with the symmetry but were not explicitly 
mentioned in Ref. [l[ . This variant of the theory is free 
of known pathologies (instabilities, overconstrained evo- 
lution, or strong coupling at low energies) which afflict 
some other variants. We will refer to this theory simply 
as Hofava gravity. 

In this paper, we consider the Hamiltonian formulation 
of Hofava gravity. While this was already considered in 
Ref. [ij (see also comments in Ref. Q), there are several 
reasons to revisit the analysis. In Ref. [||, it was argued 
that, unlike other theories with time- reparametrization 
symmetry, Hofava gravity has a nonvanishing Hamil- 
tonian. A vanishing Hamiltonian is one aspect of the 
so-called "problem of time" in quantum gravity: quan- 
tization leads to a theory in which there is no unitary 
evolution with respect to an external time Q. Here we 
find that although the Hamiltonian density is a sum of 
second-class constraints, the total Hamiltonian is a sum 
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of first-class constraints, one of which was overlooked pre- 
viously. Once this constraint is recognized, we see that 
the Hamiltonian for a closed space is indeed a sum of con- 
straints, and therefore a global version of the problem of 
time persists. 

We also consider the Hamiltonian formulation with 
asymptotically flat boundary conditions. In Ref. @ 
it was shown that in the IR limit Hofava gravity 
is closely related to Einstein-aether theory with a 
hypersurface-orthogonal aether field, in the sense that 
every hypersurface-orthogonal aether solution is also a 
Hofava solution. 1 We show that the energy of asymptot- 
ically flat solutions of Hofava gravity can be expressed as 
a surface integral at spatial infinity that agrees with the 
expression for energy in Einstein- aether theory. 



II. HORAVA GRAVITY 

In a covariant formulation [8], the dynamical objects of 
Hofava gravity are a 4 metric and a preferred foliation. In 
the Hamiltonian formulation one specifies an additional 
foliation, which need not coincide with the preferred fo- 
liation. However, it will be useful to use the preferred 
foliation in the Hamiltonian formulation, since then the 
Lagrangian depends only on first time derivatives of the 
fields. We will therefore present the theory in coordi- 
nates (t, x l ) adapted to the preferred foliation. This is the 



In fact there can be a global topological obstruction to this equiv- 
alence. If E has a nontrivial first homotopy group then it is 
possible to have an aether field that is locally hypersurface or- 
thogonal but is not normal to any global foliation. A simple 
example is provided by the "tilted aether" Bianchi type I cos- 
mologies discussed in Ref. Q, if a homogeneous aether field on 
a homogeneous spacetime with topology i? 3 X S 1 is tilted along 
the S 1 direction. 



2 



original formulation of Hofava l]. 2 In these coordinates 
Hofava gravity has the symmetry of foliation-preserving 
diffeomorphisms 



t-Hf(t), x l ^x H {x\t). 



(1) 



The metric can be expressed in Arnowitt-Deser-Misner 
(ADM) form in terms of the lapse N, shift vector N l and 
three-dimensional metric gij, 

ds 2 = N 2 dt 2 - g ij {dx i + N i dt){dx j + N 3 dt) (2) 

where letters ... = 1,...,3 denote spatial indices. 
The extrinsic cur 
and its trace K are given by 

1 



The extrinsic curvature 2£y of the constant t surfaces 



" -aJV^-Vi^-VjiVd, K = g^K ij (3) 



Ku = 



where the dot denotes partial derivative with respect to 
t, and V, is the three-dimensional covariant derivative 
compatible with the metric g^ . The acceleration of the 
congruence normal to the constant t surfaces has the spa- 
tial projection 



en = Vj In N. 



(4) 



In terms of these variables the Lagrangian density of 
Hofava gravity is 



C = 



1 



16nG H 



^gN {KijK ij - XK 2 - V{ 9ij ,<*)), (5) 



where V is a potential that depends on the 3 metric and 
the acceleration. The potential contains all terms that 
are spatial scalars of dimension up to six, where the spa- 
tial coordinates x l are taken to have dimension —1. In 
the infrared (IR) limit, the potential is dominated by the 
lowest dimension terms 



V {gn, (hi) = —£R - aa. t a l 



(6) 



where R is the Ricci scalar of the 3 metric g^ , and . . . 
refers to terms containing more than two spatial deriva- 
tives. The term V- a is also of second order in derivatives, 
but its contribution to the Lagrangian differs from that 



of a by a total derivative since iVV 



V 2 N - Na 2 



The dimensionless free parameters of the IR limit of the 
theory are a, £ and A. When a = and A = £ = 1, 
the theory reduces to general relativity. From here on we 
adopt units where 16ttGh = 1. In this paper we do not 
consider matter couplings. 



2 It can happen that such coordinates do not cover the complete 
manifold as, for example, in black hole formation by collapse [jj. 



III. HAMILTONIAN AND CONSTRAINTS 

We now consider the Hamiltonian formulation of the 
theory, following Ref. [j]. First, since the time derivatives 
of N and N l do not appear in the action, the correspond- 
ing conjugate momenta pn and pi vanish, i.e. we have 
primary constraints 

Pn = 0, Pl = 0. (7) 

The momentum conjugate to g^ is 

p« = V5(^ ~ AAV J ), (8) 
and this relationship can be inverted yielding 



K 13 = 



^0 



A 



1 -3A 



pg 



(9) 



where p = p %3 g^ . 

The Hamiltonian density H has the usual form p lJ q^j — 
C plus the primary constraints times Lagrange multipli- 
ers. Up to total derivatives it takes the form 

H = NH t + N^Hi + v i pi + vpn (10) 
Hi = -2g ik V jT P k 



ir-\-<jY) (11) 
(12) 



For now we will assume that E has no boundary, so that 
total derivative terms play no role, but they will be con- 
sidered in Sec. [V] when asymptotically flat boundary con- 
ditions are imposed. The Hamiltonian is then given by 



H 



d 3 yH(y) 



H. 



(13) 



Here and below we adopt a notation in which J means 
J s d 3 y and the dependence on spatial coordinates of in- 
tegration is suppressed. 

The primary constraints (0 must be preserved in 
time; in other words they must have a vanishing Poisson 
bracket with the Hamiltonian. Among these constraints 
there is a particular combination J NpN that generates 
constant rescalings of N. Since en is invariant under such 
rescalings, only explicit N and p^ dependence need be 
considered. Preservation of this constraint in time thus 
requires that 



(14) 



Np N ,H}= (NH t + vp N ) = 0, 



which shows that J NHt — when the constraints are 
preserved in time. It follows that the Hamiltonian is a 
sum of constraints. Moreover, the requirement that all 
constraints are preserved in time implies that the Hamil- 
tonian is first class; i.e., it has zero Poisson brackets with 
all constraints. 
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The fact that the Hamiltonian is a sum of first- 
class constraints is a feature of systems that have time- 
reparametrization symmetry. What distinguishes Hof ava 
gravity from GR is that arbitrary deformations of the 
constant-time surfaces cannot be regarded as pure gauge. 
In this regard, Hofava gravity is like a partial gauge- 
fixing of GR in which the global time-reparametrization 
freedom is left unfixed. 

Preservation of the primary constraints in time leads 
to secondary constraints 

C = SH/SN = 0, C t = 5H/5N 1 = 0, (15) 

where 

C = H t - ^V<V* (16) 

a = n h (17) 

and we have defined the vector density 

V ^ = fafe) / ^ NV - (18) 

The vector constraint Cj is the same as in GR, but the 
scalar constraint C is modified. Note that the constraint 
C is invariant under a constant rescaling of N. The 
source of this rescaling freedom is the reparametriza- 
tion symmetry (JTJ), under which the lapse transforms as 
N^N/f(t). 
In the IR limit, 

V 1 = -Za^/g^N, (19) 
and the constraint C — becomes 

There are different strategics for solving the scalar con- 
straint. For example, one could solve for the conformal 
factor of the metric while keeping the other degrees of 
freedom fixed as is often done in GR. In Hofava gravity 
there is the scalar degree of freedom TV, and it is natural 
to view C as an equation for N [U, Hj]. Note, however, 
that Eq. (|20|) can determine N at most up to a constant 
rescaling. 

The constraint equation in the IR limit can be lin- 
earized by the change of variables 3 N = n 2 , resulting in 
the equation 

Ln = (21) 
where L is the linear differential operator 

L = -4qV 2 - ~ (pV Pij + Y^P^j + (22) 



3 This was pointed out to us by M. Grillakis. 



Such an equation admits a solution if and only if the spec- 
trum of the Schrodinger-like operator L contains zero. 
Moreover, if the foliation by constant t surfaces is to be 
smooth the lapse must be positive everywhere, which re- 
quires that n(x) is positive for all x. A solution with 
positive n exists if and only if zero is the least eigenvalue 
of L: this is the familiar statement that the Schrddinger 
equation admits a unique eigenstate with no nodes, and 
this state is a ground state [l(| ■ Thus Eq. ([21) contains a 
condition on the metric and its conjugate momentum 
p 13 . If this condition is met the constraint has a unique 
(up to rescaling) positive solution for n, and hence the 
lapse is determined up to a constant scaling. 



A. Propagation of constraints 

The constraint equations C = and Ci — must hold 
at each time. One way to satisfy this requirement is to 
solve these two equations as independent constraints at 
each time. For example, one could imagine solving the 
constraint C for N as a function of the metric and its 
conjugate momentum, up to a time-dependent prefactor. 
If this could be done, the constraint would hold for all 
times and there would be no need to add an independent 
condition to ensure its preservation with time. The time 
evolution of N is governed by the Lagrange multiplier 
v, which is entirely free to begin with, so can always be 
chosen so as to produce the required time evolution of N. 
(The situation would be quite different however if V had 
no N dependence. Then the constraint would impose a 
relation between gij and p lJ alone that is not generally 
consistent with their evolution equations, except in the 
GR case.) 

Rather than solving the constraint at all times, one can 
instead solve it at one time, and then choose the Lagrange 
multiplier such that the constraint is preserved. We now 
proceed to analyze the preservation of the constraints 
in this way. This approach implements the general for- 
malism for constrained Hamiltonian systems discussed in 
Ref. [llj, and allows one to identify the first-class con- 
straints without any guesswork. 



1. Propagation of diffeomorphism constraints 

Let us consider first the diffeomorphism constraint 
Cj. Since it generates spatial diffeomorphisms, and the 
Hamiltonian is a spatial scalar (a number, not a field), 
one might think that the Poisson bracket {C,, H} would 
vanish, implying that C; is constant and therefore the 
constraint Cj = is preserved in time. There is a catch 
in this reasoning however, since C; only generates diffeo- 
morphisms of and p y , whereas the Hamiltonian also 
depends on TV and N l , so in fact {C,,i?} =^ 0. The im- 
pact of the N z dependence is transparent, since N % only 
enters H linearly in the combination N l Ci 1 and the Pois- 
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son bracket algebra of diffeomorphism generators closes, 

{a(x),C 3 (y)} = C t (y)6, 3 (x,y) + C J (x)6, t (x,y). (23) 

This N l dependence therefore produces a contribution to 
the time derivative of Ci(x) that is proportional to the 
constraint itself, which vanishes when the constraint is 
satisfied. 

The impact of the N dependence is more subtle. How- 
ever, a simple way to see that it does not spoil the con- 
servation of the diffeomorphism constraint is to modify 
the diffeomorphism constraint to include the term N t iPN 
that generates diffeomorphisms of N (and pn), which 
vanishes when the primary constraints (J7J arc satisfied 
4 . Similarly one could add the term that generates dif- 
feomorphisms of N l . The resulting extended diffeomor- 
phism constraint Ci is defined such that 

feCi = J PN C ( N + Pl C 6 N l +p ij C i g ij 



(V t N) PN +CffPi - 2g lk V,p> k (24) 



where C% is the Lie derivative. This extended constraint 
generates diffeomorphisms of all the variables, so it actu- 
ally does have a vanishing bracket with the Hamiltonian. 
Since the secondary constraints (j!5[) already imply that 
the primary constraints (O are preserved in time, the 
time independence of Ci also implies the preservation of 
the secondary constraint Ci = 0. 



2. Propagation of scalar constraint 

Next we turn to the issue of propagation in time of 
the scalar secondary constraint C = 0. First we recall 
how it works in GR, then we consider the case of non- 
projectable Hof ava gravity without the a, dependence in 
V, and finally we include the effects of this dependence. 
In GR, the potential V = —£,R is independent of N, so 
that C = Hf Moreover, these constraints are first class; 
i.e., their Poisson brackets with each other are combi- 
nations of the constraints themselves. This is obvious 
for the Poisson brackets with Ci, since Ci just generates 
diffeomorphisms. The only other bracket is [l2j 

{C(x),C(y)} = (gV^CiW+gV^CMy^y), 

(25) 

which closes on the Ci constraint (for any value of £). 
Thus the time derivative of C is a combination of the 
constraints: hence the constraints are preserved in time. 

In Hof ava gravity, let us consider first the case where 
V(gij) is independent of N. Then, when A ^ 1 and/or 



V 7^ —£R, the bracket of two C's does not close on a 
constraint, so the constraints are not first class. Preser- 
vation of C then imposes a further tertiary constraint 
that de pen ds on , p n , and N. This case was analyzed 
in Ref. [13| (see also Ref. [14|), where it was shown that 
for generic gij and p y , the only solution for N is N = 0. 
This is unacceptable since the kinetic term of the action 
§5§ is proportional to 1/N and the Hamiltonian gener- 
ates no evolution, so we will assume N ^ everywhere. 
When V = — £R, the tertiary constraint then holds if and 
only if ir/yfg is constant, and the theory is equivalent to 
GR in constant mean curvature gauge (provided such a 
gauge is accessible) [151 ] . The case where V is indepen- 
dent of N has also been considered in the presence of Re- 
type terms both in the linearized fl6l| and nonlinear fl7j 
settings, where it was found that the structure of the 
constraint algebra as well as the number of propagating 
degrees of freedom depend on which R 2 terms are in- 
cluded. In what follows we will consider the generic case 
where V depends on N and a^O. 

In the generic case, N and its spatial derivative occur 
in the ttj dependence of V(g. 



so that the constraint 
C does not commute with p^; i.e., p^ and C are second 
class. A possible way to proceed is to simply solve the 
constraint C(x) = for one of the dynamical variables, 
such as TV as suggested in Ref. [4| . The subtlety in doing 
this reduction is that iV can only be so determined up to 
an arbitrary time-dependent, spatially constant multiple. 
This is related to the fact that among the constraints 
Pn{x) and C{x) are two linear combinations that are first 
class. These first-class combinations should not be set to 
zero strongly, since the symplectic form pulled back to 
such a subspace would be degenerate, and therefore the 
Poisson brackets would be ill-defined [11]. Thus to carry 
out the reduction we identify the full set of first-class 
constraints. 



B. First-class constraints and Hamiltonian 

For convenience in the analysis, we write the Lagrange 
multiplier v in the form 



v = Nw + N l V t N, 



(26) 



where the function w is initially arbitrary. In terms of w 
the Hamiltonian density (TTU1) takes the form 



H = NH t + N l d + v l pi + Nwp 



N 



(27) 



4 Such a modified constraint was considered in Ref. Q but not 
explicitly identified as a generator of diffeomorphisms. 



up to a total derivative. The condition that C be pre- 
served in time is then 

{C(x), H] = J N({C(x), H t } + w{C(x), PN }^ = 

where we have dropped the term proportional to {C l , C}, 
which vanishes when C = 0. By solving for w, we ensure 
that the constraint C = is preserved in time, assuming 
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it has been solved at some initial time. Thus rather than 
solving the nonlinear equation C = for N separately 
at each instant of time, one can instead solve a linear 
equation for w. The presence of the free function w can 
solve the overconstraining problem, provided the bracket 
(|28| can be set to zero by solving for w. 

The Lagrange multiplier w is not completely deter- 
mined by preservation of C, instead the Eq. (|28|) deter- 
mines w up to a solution of the homogeneous equation 



Nw{C(x), PN } = 0. 



(29) 



Each such solution is a gauge freedom in the evolu- 
tion and corresponds to a primary first-class constraint, 
J wNpN- In this case w{x) = k is such a solution, where 
k is constant. The existence of this solution follows from 
the fact that the constraint 



Np 



N 



(30) 



generates constant rescalings of N, and C is invariant 
under such a rescaling. 

Equation (|28[) is a linear partial differential equation 
for w of up to sixth order. In the IR limit it reduces to 
a second-order equation, 



(31) 



This equation is elliptic, and in fact the left-hand side can 
be written as y/gW 2 w where V is the covariant derivative 
compatible with the metric gij — N 2 gij. It follows that 
the solution exists for w provided the integral of the right- 
hand side is zero, 



NC, / NU t \ = 



(32) 



which follows immediately from the definition of C (TH 
This solution is unique up to the addition of the homo- 
geneous solution w{x) = k. 

Beyond the IR limit, the equation determining w is of 
higher order and might admit nonconstant homogeneous 
solutions. However any theory admitting additional so- 
lutions to Eq. (|2"9"|) will have additional first-class con- 
straints, which may be regarded as generators of gauge 
transformations 11] . For a generic choice of parameters, 
we expect there will be no additional gauge symmetry, 
so that the general solution of (|2"51) has the form 

(33) 



w = w[g 2J ,p lJ ,N] + k, 



where k is constant and w is a particular solution of 
Eq. (JMl). 

The solution for w can be substituted back into the 
Hamiltonian to obtain a Hamiltonian that preserves the 
constraints in time, 



H = J N(H t + wp N ) + K J 



t + wp N ) + k / N PN + / N l C % + v % p. 



By varying the Hamiltonian by a constant rescaling of 
N, we can see that the first two terms in Eq. (jM)) are the 
Poisson brackets of the Hamiltonian with the first-class 
constraint f N P n, 



H, / Np 



N, 



N{U t + wp N ). 



(35) 



The Hamiltonian H is therefore a sum of first-class con- 
straints generating the two types of foliation-preserving 
diffcomorphisms: global time reparametrizations, and 
spatial diffeomorphisms. 

Observables in Hof ava gravity must be gauge-invariant 
functions; that is they must commute with all first-class 
constraints. Since the Hamiltonian is a sum of first- 
class constraints, observables must also have zero Poisson 
bracket with the Hamiltonian and hence be conserved in 
time. For example, the volume of the spatial slice E 
is diffeomorphism-invariant, but is not time-independent 
since it fails to commute with the constraint J NC: 



{J 14-^1 



Np. 



(36) 



(34) 



Hence the volume is not an observable in Hof ava gravity, 
contrary to what was claimed in Ref. [H . This reflects the 
fact that it is meaningless to label an observable by the 
"f coordinate in a theory that has t-reparametrization 
symmetry. 



IV. REDUCED PHASE SPACE 

We have expressed Hofava gravity as a Hamiltonian 
system with second-class constraints. In the presence 
of second-class constraints it is possible to reduce the 
phase space by solving the constraints for one or more 
of the dynamical variables as a function of the others. 
The Hamiltonian and symplectic form are then restricted 
to the reduced phase space, and the restriction of the 
symplectic form defines a nondegenerate Poisson bracket 
on the reduced phase space known as the "Dirac bracket" 

In Refs. @, 3, it was proposed that the constraint 
C = (fT5j) be solved for N. There are two related is- 
sues with this strategy. First, N can be determined only 
up to a time-dependent rescaling. Second, among the 
constraints are two first-class constraints; setting these 
to zero strongly would result in a degenerate symplectic 
form and therefore an undefined Dirac bracket. 

To reduce the system, we need to impose as many 
linear combinations of the constraints as possible with- 
out setting the Hamiltonian to zero. In order to have 
a set of constraints that determines N completely, we 
choose the gauge in which the average lapse is one. Once 
this new constraint is introduced the only first-class con- 
straint that remains (apart from the diffeomorphism con- 
straints) is the Hamiltonian. We therefore impose the 
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following constraints strongly, 



/ N^g = J y/9, Pn = 0, C = C 0y /g 



(37) 



and eliminate N from the phase space. Here Co is a 
function of time whose presence is necessary to ensure 
that the Hamiltonian is not set to zero strongly. Co can 
be expressed in terms of N by integrating the identity 
C = Coy/g, or alternatively, since TV is completely deter- 
mined by the constraints, Co can be written in terms of 
the metric variables gij , p %3 . 

The reduced phase space has coordinates gij,p lJ and 
its dynamics are expressed in terms of the Hamiltonian 
and the Dirac brackets [llj. Although generically the 
Dirac brackets are different from the Poisson brackets, in 
the case where the second-class constraints can be writ- 
ten in the form N — f[gij,p % ^\, Pn = 0, the Dirac and 
Poisson brackets between gij and p 13 coincide. Since the 
constraints (|37|) do not restrict the values of g^ and p %3 , 
they are of this special form and the Dirac bracket {•,•}* 
on the reduced phase space is the same as the canonical 
Poisson bracket, 

{9i^),P kl {y)Y = {g, 3 {x),p kl {y)} = S\ k 5 l j5{x,y). (38) 

The Hamiltonian on the reduced phase space is found 
by substituting the second-class constraints (|37| into the 
Hamiltonian (1341): 



(39) 



where V is the volume of the spatial slice E. Note that 
the terms involving w are absent, since their only role was 
to preserve the constraint C = 0, which always holds on 
the reduced phase space. The equations of motion are 
then 



H = VC a + / (N'd+v'p. 



V 



SCo 
Sp 13 
SCo 



-V 



Sgi 



(40) 
(41) 



where we have used the fact that Co = on shell. 

In the IR limit, —Co can be expressed as the smallest 
eigenvalue of a linear operator, 



Ln 



-C n, 



(42) 



where L is given by Eq. (|22j) . and n = v N. The time evo- 
lution then depends on first-order changes of the eigen- 
value Co with respect to g, and p. These can be found 
by first-order perturbation theory, e.g. 



„SC S _ 
V- — = —z — / y/gnLn 
ogij ogij 



(43) 



and similarly for p. This formula can then be rewritten 



in terms of N 
.SCo 



2N 



Sp i3 ~ 
,SC_o 
Sgij 



x 



pgtj 



v 

Vp^ = y/gN (aa k a k + £R) g 



1 - 3A 



(44) 



N f ki 



X 



3A 



p 2 ) g 13 



aiy/gNa 1 



2N 



P' k Pk 



PP- 



yfg V rK ' 1 - 3A 

+ i^fg (V*V J - R ij - g ij V 2 ) N. (45) 

These equations are equivalent to those obtained by vary- 
ing the Lagrangian ([5]). 



V. ASYMPTOTICALLY FLAT CASE 

In the preceding discussion, it was assumed that the 
spatial manifold E is compact with no boundary. When 
there is a boundary (or in the case of asymptotic flatness, 
an asymptotic region with prescribed falloff conditions), 
variation of the Hamiltonian will result in a total deriva- 
tive that can be written as a boundary term. To have a 
well-defined variational principle, additional terms must 
be added to the Hamiltonian to make this boundary term 
vanish. The appropriate boundary terms in the asymp- 
totically flat case are determined by the falloff conditions 
on the fields, 

9ij = Sij+Oir- 1 ), Ka = 0(r- 2 ), R = 0(r~ 3 ) (46) 

where r is the radial coordinate of a coordinate system 
in which the metric asymptotically approaches the Eu- 
clidean metric 5y . In order for the Hamiltonian to define 
an asymptotic time translation, we must have N l — > as 
r^oo. and we choose a falloff on N such that 



N=l + 0(r- 1 ), ViN = 0{r- 2 ) 



(47) 



as r — > oo. The contributions to the boundary term Hg 
in the Hamiltonian then come only from the potential 
V. For any term in V that is fourth order or higher 
in derivatives, the corresponding boundary term is of at 
least third order in derivatives and hence vanishes subject 
to (I46p . The remaining terms in V, which are of second 
order in derivatives, are R and a 2 . The a 2 term does not 
contribute a boundary term because aiSN = 0(r -3 ), and 
the variation of R with respect to g^ can be canceled by 
adding to the Hamiltonian the boundary term 



Hd = Zf {dig ]k - d jgik )S ik n 3 



(48) 



d 2 ay/h where a are coordinates on the 
sphere at infinity 9E, h is the induced metric on 9E, and 



Here / - :r . is _. 
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n l is the outward unit normal. This term is nothing but 
the usual boundary term for general relativity, the ADM 
mass, with a factor of £ coming from the action ([5]). We 
have therefore obtained a Hamiltonian that leads to a 
consistent variational principle in the asymptotically flat 
setting, given by 



H= H + H dl 



(49) 



where the first term is given by (| IQtil 3|) . 

The primary and secondary constraints (|7|) and (|15j) 
remain the same in the asymptotically flat setting. The 
difference comes when solving for the undetermined part 
of the Lagrange multiplier w to ensure the propagation 
of the constraint C (|29|) . In the asymptotically flat set- 
ting there is no global first-class constraint / Np^, when 
the global time-reparametrization symmetry is broken by 
the asymptotic value of N. Formally, the first-class con- 
straint is absent because in order for the flow generated 
by J wNppj to preserve the boundary condition (|4T|) . wc 
must have w(x) — > as x — > oo. Therefore the homoge- 
neous solution w = 1 to Eq. (f29|) is not admissible and 
the constraints pm and C are purely second-class. 

The Hamiltonian can be expressed as a sum of bulk 
constraints f| 13[) and a boundary term by multiplying the 
definition of C (IT5|) by N and integrating 



(50) 



Terms in V with four or more derivatives contribute 
terms in V 1 with three or more derivatives, and these 
vanish subject to the falloff conditions (|4"6")l so that we 
can replace V % with its IR limit, Eq. <jTUJ) . Combining 
these results with (|2"T)) we find that the Hamiltonian is 



NU t = NC+ & n t V\ 



H 



(iVC + wN PN + N l C l + v*p^ + S, (51) 



which is a sum of local constraints, plus a boundary term 
£ , the total energy. The total energy £ is the on-shell 
value of the Hamiltonian, which is the sum of two bound- 
ary terms 



£ = £, <j> {dig-jk - ihn.i. •.•>'"'•//•' - 2a I n l d l 



N. 



(52) 



We emphasize that this is the total energy for the full 
Hofava gravity, not just in the IR limit. 5 

In the limit of large r, the metric can be treated as a 
perturbation of flat spacetime. In Newtonian gauge, 



N = 1 + V, 9ij = (1 - 24>)8. 



(53) 



5 An expression for the total energy in Horava gravity was recently 
found in Ref. [19|| using the Noether charge formalism. 



the energy takes the form 

£ = (bn i di(^-2onp). 



(54) 



Note that the ADM energy in GR depends only on <\> 
whereas in Hofava gravity there is also dependence on ij) 

In the static spherically symmetric solutions, 4> = 
and the asymptotic form of the solution to 0(1 /r) is 



N = 1 



2r' 



9ij= (l + y^Sij, 



(55) 



where r$ is a parameter with dimensions of length. Re- 
turning to units in which IQ-kGh ^ 1, the energy of this 
solution is 



£ 



2G [ ^ 2 



(56) 



The gravitational constant that appears in the Newto- 
nian limit is not Gh but Gn = Gh/{£, — \oi) [3, 
Thus, if we identify tq = 2G^M in the weak field limit 
then £ = M. 

Horava gravity in the IR limit can be thought of 
as the hypersurface-orthogonal restriction of Einstein- 
aether theory, where the restriction is at the level of the 
Lagrangianp, [H Because of this close relation, it 

is not surprising that the total energy expression (|54[) is 
identical to the expression for total energy in Einstein- 
aether theory found in Ref. .21], when V^lniV is iden- 
tified with the acceleration of the aether, and when the 
parameters of Hofava gravity are identified with the pa- 
rameters of Einstein-aether theory as in Ref. @ . 

The reduction of the phase space with asymptotically 
flat boundary conditions proceeds in much the same way 
as in the compact case. The difference is that the con- 
straints C are all second class, so the equation C = 
can be solved strongly as a differential equation for N 
with the boundary condition N — > 1 at spatial infinity. 
The reduced phase space has coordinates gij , p lJ with the 
canonical Poisson brackets (I3"51) . The Hamiltonian on the 
reduced phase space is 



H 



NC+ / (N i C i + v i p i )+£, 



(57) 



where N is treated as a functional of gij and p 13 defined 
by the constraint. 



VI. CONCLUSION 

We have derived the Hamiltonian formulation of 
Hofava gravity for both closed and asymptotically flat 
spatial manifolds, extending the analysis initiated in 
Ref. Q. In contrast with GR, the scalar constraints are 
not generators of surface deformations. These constraints 
are of second class, except for a single linear combination 
that is first class, which generates transformations be- 
tween leaves of the preferred foliation. The Hamiltonian 



8 



is a sum of constraints that generate three-dimensional 
diffcomorphisms and global reparametrizations of time. 

We have also considered a phase space reduction of 
Hofava gravity in which the second-class constraints are 
solved formally for the lapse. This is complicated by 
the fact that the constraints have a global first-class part 
that cannot be set to zero strongly. The second-class con- 
straints take the form of a nonlinear differential equation 
for N, but in the IR limit they are linear in \/~N. In the 
IR limit, the Hamiltonian is expressible in terms of the 
least eigenvalue of a certain differential operator, and is 
constrained to vanish. 

In the asymptotically flat setting, the timc- 
reparametrization symmetry is fixed by the constant 
value of the lapse at infinity, so the associated global 
constraint is absent. Instead the Hamiltonian acquires 
additional boundary terms that contribute to the 
total energy. The energy depends on both Newtonian 
potentials <f> and ip, not just on </> as in GR, in exact 
agreement with the result from Einstein-aether the- 



ory. For the static spherically symmetric solutions, 
4> = %jj and the energy takes the same form as GR, 
with a rescaled Newton constant. We do not know 
whether there are solutions in which <fi ^ ip, for which 
the energy would not be just a rescaling of the GR energy. 

Note added: Shortly after version one of this pa- 
per appeared on the arXiv, J. Bellorin and A. Restuccia 
posted Ref. [22j], which treats many of the same issues 
addressed here. 
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